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We generalize an analogy between rotating and stratified shear flows. This analogy is summarized 
in Table 1. We use this analogy in the unstable case (centrifugally unstable flow v.s. convection) 
to compute the torque in Taylor-Couette conflguration, as a function of the Reynolds number. At 
low Reynolds numbers, when most of the dissipation comes from the mean flow, we predict that the 
non-dimensional torque G = T/v^L, where L is the cylinder length, scales with Reynolds number 
R and gap width rj, G — 1.4677^''^(1 — rj)~'^^'^R?^'^ . At larger Reynolds number, velocity fluctuations 
become non-negligible in the dissipation. In these regimes, there is no exact power law dependence 
the torque versus Reynolds. Instead, we obtain logarithmic corrections to the classical ultra-hard 
(exponent 2) regimes: 

^ = 50^^^ 

(1 - 7?)3/2 ln[r;2(l - 77)i?Vl04]3/2 • 

These predictions are found to be in excellent agreement with available experimental data. Predic- 
tions for scaling of velocity fluctuations are also provided. 

PACS numbers: 47.27 -i Turbulent flows, convection and heat transfer - 47.27.Eq Turbulence simulation and 
modeling - 47.27.Te Convection and heat transfer 



I. MOTIVATION AND OBJECTIVES 

At sufRciently large Reynolds number, the fluid be- 
tween co-rotating coaxial cylinders becomes turbu- 
lent, and a significant momentum transport occurs 
between the two cylinders. In the case with ro- 
tating inner cylinder and resting outer one (the so- 
called Taylor-Couette flow), detailed measurements 
show that the torque applied at cylinders by the tur- 
bulent flow is a function of the Reynolds number R. 
There is no clear consensus about this dependence 
yet: marginal stability computation of King et al 
or Barcilon and Brindlcy [11^ predict that the non- 
dimensional torque G = Tjv^L, where L is the cylin- 
der height should vary like G-i?^/3. Old experimen- 
tal data indicated the existence of two scaling regimes, 
one for R > 10'' where the exponent is 1.5, and one for 
larger Reynolds number, where the exponent switches 
towards 1.7— 1.8 [T51 [301 IS] • Recent high precision 
experimental data yielded no region of constant expo- 
nent, and revealed a transition with a marked change 
of approximate slope of G as a function of i? [SI [13] . 
This observation led Eckhardt et al [T^ to propose 
a new theory, in which the dependence G versus R 
is through a superposition of scaling laws (describing 
contribution from a boundary layer and the bulk flow) . 
They claim that this superposition flts the data better 
than the Prandtl-Karman skin friction law proposed 
by [Sj HI [13]. Note that all the scalings are within 
the theoretical bound derived by Doering and Con- 
stantin [T], which implies that the non-dimensional 
torque cannot increase faster than R^. 



The observational features are reminiscent of heat 
transport in turbulent thermal convection, where ap- 
proximate scaling laws and transition between differ- 
ent regimes have also been observed (for a review see 
[TO]). In fact, this similarity is pointed out in j^, and 
in [T7] , and similar techniques are used in [TH j^H] and 
in [ini [17] to derive theoretically the scaling regimes 
in the Rayleigh-Benard system, and in the Taylor- 
Couette system. However, the similarity is more than 
superficial: as well known since several decades [3], 
there is an exact analogy between equations of mo- 
tions of rotating and stratified shear flows (stable or 
not). There must therefore exist an exact analogy be- 
tween the momentum and heat transport in these two 
systems, although it has so far never been explored. 
Our goal here is to derive this analogy, and examine 
its consequences in the unstable regime, where the an- 
gular momentum or temperature stratification leads 
to a linear instability. We thus in this paper mainly 
focus on the analogy between centrifugally unstable 
Taylor-Couette flow, and convection. 



II. THE ANALOGY 

A. Reminder 

The root of the analogy can be found in the Lamb 
formulation of the Navier-Stokes equations: 

9tu - u X cj = -V (^p + ^) + j^Au, (1) 
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where w is the vorticity, p is the pressure and v the 
molecular viscosity. The constant density has been set 
equal to one for simplicity. In a rotating shear flow 
u — V{r)eg, the vorticity is only in the axial direction 
and the Lamb vector u x w acts only in the radial 
direction. Its contribution can be split in two parts: 
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Gr • (U X = VS + 2-C, 
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where C = rV is the angular momentum, f2 = V/r 
is the angular velocity and S = rdr^ is the shear. 
The first contribution is the exact analog of the Lamb 
vector of a pure shear flow, in a plane parallel geom- 
etry. The second contribution reflects the stabilizing 
influence of the Coriolis force. Its analog would be 
produced by temperature stratification in the span- 
wise direction of a planar shear flow. The equation 
of angular momentum conservation then suggests to 
split further this analogy by requiring that 2fl/r ~ j3g, 
where g is the gravity, and (3 a coefficient of thermal 
expansion, and C ^ 0, the potential temperature. 
This remark is at the heart of the analogy between 
the stability properties of rotating and stratified shear 
flows, and has been used in the past (see e.g. 
Our point here is to show that it can be extended into 
the turbulent regime, via a new Langevin model of 
small-scale turbulence. This new model is based on 
Rapid Distorsion Theory, i.e. on linearized equations 
for the small-scale motions. This linear structure ex- 
plains the possibility of extension of the analogy to- 
wards the turbulent regime. 



B. The turbulent model 

The turbulent model has been described and tested 
in for general 3D flows, in [351 US for shear flows 
and in [2Tti24] for stratifled shear flows. In this model, 
the dynamics of the turbulent flow is obtained from 
solutions of two coupled sets of equations. The flrst 
one described the dynamics of the mean velocity U: 



dtU, + d,U,U,+d, <uW. > 
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-d,P + i^djd,U,, (3) 



Here, the primes denote fluctuating quantities and <> 
the averaging. To close the system, we need < u'^u'j >. 
They are obtained as solution of a linear, stochastic 
equation valid for localized wave-packets of velocity 
and temperature: 



Dtiii 



0, 



(4) 



where 



M(x,k,t)= y"g(|x-x'|)e*''-(^-^')u(x',t)dx', (5) 

g being a function which decreases rapidly at inflnity. 
We have dropped primes on fluctuating quantities for 



convenient notations and introduced the total deriva- 
tive Dt = dt -\- Ujdj — dj{Uiki)dkj ■ Once the solu- 
tions of Q have been computed, the Reynolds stress 
is found by an inverse Gabor transform as: 



<u[u'.j>= J dk (u-(k,x,t)^(-k,x,t) -I- (u-(-k,x,i)?^^(k,x, i)) 

(6) 

Note that the linear part of Q is exact and describes 
non-local interactions between the mean and the fluc- 
tuating part. The major approximation of the model 
is to lump the non-linear terms describing local inter- 
actions between fluctuations into a turbulent viscosity 
i^f. The force / appearing in Q is a small scale ran- 
dom forces which is introduced to model the seeding 
of small scales by energy cascades (for example via 
turbulent structures, detaching from the wall). 



C. The Taylor-Couette case 

In the Taylor-Couette (rotating shear flow) case, the 
equations for the azimuthal component of the velocity 
V{r) simplify into: 



dtV + <uv>^v (v'^V - ^ 



(7) 



The equation for the fluctuations {u,v,w) become: 
Dtu = 2^ {n + S)u + 2nv (l-^)- + 
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DtV = 2'^u{n + S)-2'^vn-{2n + S)u-Vtk^v + fe, 
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Dtw = 2^^u{Q. + S)-2-^vVl-Vtk^w + f^. 



fc2 



(8) 



Here, we have used the incompressibility to eliminate 
the pressure. These equations have to be supple- 
mented by the equations describing the ray trajec- 
tories: 

f = 0, 9 = 0., i = 0, 
fcV = -kgS, ke = 0, 4 = 0. (9) 

We now introduce a pseudo-temperature 

6 = r{v — -j^w) — ir-r^, (10) 

where is the radial vorticity. With this temperature 
and using the incompressibility condition k-u = 0, we 
can rewrite (jsl) as: 



Dtu = 2S^u + 2-^e^vteu + U 
k^ r k^ 

Dt0 = {2n + S)ru~Vtk^e + f0, 



(11) 
(12) 
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The set of equation (12 1 is the exact analog of the 
equations describing the behavior of vertical and tem- 
perature fluctuations in a stratified shear flow (see 
[3TH23] for their expression), provided the correspon- 
dence summarized in Table 1. holds. Note that the 
analog of the temperature is not the angular momen- 
tum, but related to the z-integral of the radial vor- 
ticity (in Gabor variable, integration on z is done via 
division by kz)- At large scale, since the velocity pro- 
file is axi-symmetric, this integral of the radial vortic- 
ity reduces to the angular momentum, as previously 
suggested [1]. 



D. Stability and Importance of axi-symmetric 
modes 



The correspondence described in Table 1. gener- 
alizes the well-known analogy established previously 
[4] for the stability analysis under axi-symmetric per- 
turbations (case where kg = 0). In particular, from 
( [T2| ), one can write a compact differential equation for 
u! = uk'^ : 



Dfuj 
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fc2 



Vk^i^ + = 0, 



(13) 



where = 2^{2Vl + S) is the epicyclic frequency. Us- 
ing the ray equation ([9]), one can then define a non- 
dimensional number B = K^k'^/{kg + k^) (the Brad- 
shaw number) which governs the stability of the wave 
packet along the trajectory. This number is the ana- 
log of Richardson number in stratified shear flow. For 
example, it can be shown that in the absence of diffu- 
sion, the amplitude of the wavepacket has a monotonic 
(growing for one mode, decaying for another one) be- 
havior at late time for B < 1/4, while it becomes os- 
cillatory for i? > 1/4. Clearly, the oscillatory behav- 
ior creates dephasing effects for the Reynolds stresses, 
which may lead to its pure cancellation, thereby re- 
moving the influence of the small scales onto the large 
scale. We therefore identify the regime with B > 1/4 
as a regime with purely laminar motions, where tur- 
bulence effects are strongly suppressed. This property 
tends to favor bi-dimensional modes (those for which 
kz = 0) since in this case the epicyclic frequency can 
take any value for non-oscillatory behavior. The in- 
clusion of diffusion changes the mode selection. One 
can indeed check that the viscous decay is propor- 
tional to exp{—R^^tS) rather than exp(— i?^^(i5)^) 
for non-axi-symmetric perturbations. This shows that 
axi-symmetric perturbation (with kg = 0) are favored 
with respect to non-axi-symmetric perturbation. In 
this case, the Bradshaw number becomes independent 
of the wavenumber of the wave packet, and one can 
identify a new boundary of stability according to its 
sign: when it is positive, axi-symmetric perturbations 
can be exponentially amplified and superseed viscous 
decay. We call this regime "unstable". It is the ana- 
log of the convective regime in the stratified case. In 



the sequel (Section 3), we shall concentrate on this 
regime, leaving the other regime for further study. 



E. Completion of the analogy in the unstable 

case 

For axi-symmetric modes, uj^/kz — v and the equa- 
tion for the mean angular momentum C can then be 
written in equivalent form: 



1 / C 

dtC + -dr(r <u0>)=iy{ W^£/r - 



(14) 



Comparing this equation with the equation giving the 
mean temperature profile, we finally remark that the 
only difference lies in the viscous terms, because in 
cylindrical coordinates, the Laplacian includes terms 
describing curvature effects. In the most general case, 
this forbids the analogy to be drawn at the level on 
mean profile (ie after integration over r ofeq. 0: 
for example, it is well known that in stratified shear 
flow, the laminar temperature profile is linear, while 
its analog, the laminar angular momentum profiles 
varies like: L ~ Ar^ + B. In many Taylor-Couette 
experiments, however, the gap between the two cylin- 
ders is small, and curvature effects can be neglected. 
One can for example check that the angular momen- 
tum in the experiments by [S] is linear in the laminar 
regime, while it fiattens at the center of the gap in the 
turbulent regime, exactly like its temperature analog. 
In the sequel, we shall assume a small gap geometry, 
and neglect curvature effects. 



III. APPLICATION OF THE ANALOGY IN 
THE UNSTABLE CASE 

The present analogy is the turbulent generaliza- 
tion of a previously known analogy for axi-symmetric 
modes. In the sequel, we shall use previous consid- 
erations about stability of axi-symmetric modes to 
assume that the turbulent properties are dominated 
by the contribution of the axi-symmetric modes, i.e. 
restrict ourselves to these modes. The relevence of 
this approximation will be tested by comparisons of 
its predictions regarding some characteristic quanti- 
ties measured, in the turbulent regime with experi- 
mental data. 



A. Stability 

In the unstable regime, a classical parameter de- 
scribing the intensity of the convection is the Rayleigh 
number: 



Ra 



(15) 
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where D is the size of the ceh in the stratified direc- 
tion, AO is the temperature gradient apphed. In most 
convection experiments, this number is unambigu- 
ously defined because of the constancy of the gravity 
at the scale of the experiment. In the Taylor-Couette 
case, the gravity depends on the perturbation, and 
one may wonder how to define this Rayleigh number 
in a general way. In a recent analysis of stability of 
Taylor-Couette experiment, Esser and Grossman [3] 
suggested to evaluate this factor at the gap center, 
''c = ('^1 + '^2)72, leading to 



Z Orl^ o \r—rc ; 
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(1 — ?7^)^ r. 



(16) 



In the sequel, we shall use a star label to refer to 
analog quantities. In (16), we have used = v 
and introduced the Reynolds number R = ridili/v, 
where ri is the internal radius, ili the rotation rate 
at the inner radius, d the gap width and 77 = rxjr^. 
Note that the analog Rayleigh number i?a* varies with 
the radial aspect ratio 77. In the small gap approx- 
imation 77 — > 1, experiments show that the critical 
Rayleigh number tends to a constant i?a* « 1706. 
This value is very close to the value Ra = 1707.762 ob- 
tained in Rayleigh-Benard convection for rigid bound- 
ary conditions [5]. In sheared convection, the criti- 
cal Rayleigh number is modified with respect to this 
theoretical value, and display corrections quadratic 
in the Reynolds number based on the shear. In 
the present case, these correction are proportional to 
(1 — riY 0, and the critical Rayleigh number stays 
close to the un-sheared value Ra — 1707.762. A last 
modification of the critical Rayleigh number occurs 
because of lateral wall effects. As a result, the critical 
Rayleigh number increases with decreasing aspect ra- 
tio r (lateral width over radial width). For example, 
for r = 5,2,1,0.5, Ra^ = 1779,2013,2585,12113. In 
most Taylor-Couette experiments, the aspect ratio is 
very large (typically above 8 or so) . So the analog crit- 
ical Rayleigh number is close to 1708. However, many 
modern convection experiment (reaching very large 
Rayleigh numbers) deal with a rather small aspect 
ratio (r ^ 1). This unfortunately limits the possibili- 
ties of direct comparisons between the Taylor-Couette 
experiments and the convective experiments to values 
close to the onset of instability. For larger values of 
Rayleigh numbers, we shall use extrapolations. 



B. Angular momentum transfer 

A second interesting quantity in convection is the 
non-dimensional heat transfer Nu — Hd/KAT, where 



H is the heat transfer. Via the analogy, the analog of 
this is the non-dimensional angular momentum trans- 
fer, which can be computed using the non-dimensional 
torque G = T/v^L, where L is the cylinder length: 



G 
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G (l+^)(l-^)2 

R Aurj 



(17) 



The normalization by Giaminar ensures that in the 
laminar case, Nu^, = 1, like in the convective analog. 



1. Instability onset 

Theoretical [7] and experimental [S] studies of con- 
vection near threshold lead to identification of two 
regimes just above the critical Rayleigh number: 
for e = Ra-Ra^ < 1, a linear regime in which 
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Ra 
Ra, 



(18) 



The constant Ki depends on the Prandtl number. 
Pr = 1, it is Ki « 1/0.7 = 1.43 [7j. 

for larger e, a scaling regime in which j6] 



For 



{Nu - 1) 



Ra 

Rac 



^1-23 



(19) 



Here, K2 is a constant which is not predicted by the 
theory. In Fig. [l] we show how the results of Wendt 
obtained with 77 = 0.935 near the instability threshold 
compare with these two predictions. One sees that the 
linear regime is indeed obtained for e < 10, while the 
scaling regime is obtained for larger values of 10 < 
e < 100. Further from the threshold, one needs to 
compare with the turbulent theories of convection. 



C. The classical turbulent regimes 



Using (17 1, we can single out some interesting 
regimes. Note that since k* = i^, we are in the case 
of unit Prandtl number convection, i.e. for exam- 
ple convection in Helium {Pr > 0.7). In the classi- 
cal theory of convection, one usually considers three 
regimes: a first one, labeled as "soft turbulence", in 
which Nu - i?ai/3, 5 X 10^ < i?a < 2 X 10^ [H]; then 
for 2 X 10^ < Ra < 10^^, a "hard turbulence" regimes 
in which Nu ^ Ra^/"^ IT; finally for Ra > 10^\ a 
ultra-hard turbulence regime in which Nu ~ Ra^^^ 
[15, 27 . Using the analogy, we see that these three 
regimes translate into: for 707 < Rri{l-ri)^/'^ < 4472, 
G ~ i?^/^ for 4472 < i?77(l - tjY^^ < 3 x 10^ 
G - i?"/^; for (1 - riy/^r]R > 3 x 10^ G - R^. 
To evaluate the boundary between the two regimes, 
we have used (16) at 77 = 1. 
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1 10 

E = (Rfl-Ra,)/Ra, 



FIG. 1: Comparison of the theoretical near instability 
onset behavior with the data of Wendt [181. The symbols 
are the experimental measurements. The two lines are the 
theoretical formula predicted by analogy with convection 
for e < 1 {{Nu - l)Ra/Rac = 1.43e and for e > 1 {{Nu - 
l)Ra/Rac ~ e^'^^.) In the latter case, the proportionality 
constant is not constrained by the analogy, and needs to 
be adjusted for a best fit. 



The first regime has been predicted by [TTl [12] us- 
ing marginal stability analysis. The third regime can 
be derived from Kolmogorov type arguments (see e.g. 
[5]). It also corresponds to some upper-bound in the 
angular momentum transport [1]. The intermediate 
regime is new, and leads to a scaling exponent of 
1.57. Experimentally, some of this scaling regimes 
have been reported, but not in the same sequence: in 
his experiments with 0.680 < 77 < 0.935, Wendt [E] 
reports a scaling exponent of 1.5 for 400 < R < 10^, 
followed by a scaling exponent 1.7 for 10** < R < 10^. 
In more recent experiments, Lathrop et al [S] mea- 
sure a "local" exponents dln{G)/dln{R) which varies 
continuously from 1.2 to 1.9, with a transition at 
i? ~ 1.3 X 10'' (for 77 = 0.7246). This transition was 
later found to correspond to a modification of coher- 
ent structures in the fiow 1 1 3_j . Remarkably enough, 
the analog Rayleigh number characterizing this tran- 
sition is Ra^, = 2 X 10^, like in convection. We may 
therefore interpret it as the boundary between "soft" 
and "hard" turbulence. However, the scaling reported 
in [5] does not seem to fully correspond to the soft and 
hard turbulence scaling. In the sequel, we wish there- 
fore to explore a new possibility, based on logarithmic 
corrections to scaling. 



tion [TB]. In a recent work, we used the turbu- 
lent model to analytically compute the heat trans- 
port in a convective cell. At Pr = 1, we found 3 
different regimes: at low Rayleigh number, the dis- 
sipation is dominated by the mean flow, and Nu = 
KiRa}/'^Pr^^/^^\ at larger Rayleigh number, the ki- 
netic energy dissipation starts being dominated by ve- 
locity fluctuations, and the heat transport becomes 
NuPr^l^ = X2i?ai/^/ln(i?aPr2/3/20)2/3. Finally 
at very large Rayleigh number, the heat dissipation 
becomes also dominated by (heat) fluctuations, and 
Nu = KsRa^/^/lniRa/Racf/^. Fig. [2] shows the il- 
lustration of these 3 regimes in a Helium experiment 
of [27], with the three fits corresponding to these 3 
regimes. From this graph, we obtain Ki = 0.31, 
K2 = 0.45, K3 = 0.023 and Rgc = 2 x lO'^, for an 
aspect ratio of 0.5. These constants tend to decrease 
slightly for larger aspect ratio by an asymptotic fac- 
tor of about 0.75 (at Ra = 10*^, see table 1 of [17]). 
The small aspect ratio of the experiment also increases 
the critical Rayleigh for instability from near 1708 to 
near 4 x 10**. The boundary between the regime 1 
and 2 lies at Ra = 1.5 x 10^. It is characterized 
by a change in the temperature statistics, going from 
nearly Gaussian to exponential. The boundary be- 
tween the regime 2 and 3 is somehow ill defined, and 
lies between i?a = 2 x 10^° and Ra = 10". Note 
than in a similar experiment, ran by another group, 
the third regime was not detected, even at Ra = 10^^ 
[5D]. The reason of this difference is not yet known. 
A possibility would be that different boundary condi- 
tions may or may not allow the growth of the temper- 
ature perturbation, thereby favoring or inhibiting this 
last regime [25] . 

The translation of the three logarithmic regimes us- 
ing the analogy gives a priori three possible regimes 
in the Taylor-Couette experiments. 

In the regime 1, we get: 
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D. The logarithmic turbulent regimes 

The observation by [8 that no scaling prevails 
for the angular momentum transport has in fact 
its counter-part for the heat transport in convec- 



In these expressions, we have introduced 5 unknown 
coefficients, which a priori depend on the aspect ra- 
tio. Since there is no available large Rayleigh number 
large aspect ratio convection experiments, we shall ex- 
trapolate or fit these coefficients by comparison with 
Taylor-Couette data. 
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FIG. 2; Illustration of the three scaling regimes found 
in convection in Helium for Nusselt vs Rayleigh. The 
symbols are experimental measurements by [27]. The 
lines are theoretical prediction by using an analyt- 
ical model of turbulent convection. "Soft" turbulence 
regime(mean flow dominated): power law Nu ~ Ra^^* 
(full line); "Hard" turbulence regime: (velocity fluctua- 
tion dominated) Nu ~ Ra^^^ / {ln{Ra)f^^ (dotted hue); 
"Ultra-hard" turbulent regime: (temperature fluctuations 
dominated) Nu ~ i?a^''^/(ln(i?a))^/^ (dashed line) 



E. Comparison with experiments 

For this, we use torque measurements from Wendt 
[H] and [HI [T^. The regime 1 should be observed 
at rather moderate Reynolds numbers. Therefore, it 
explains very well the old measurements by Wendt 
[18] who found the same exact dependence in rj and 
R for 400 < R < 10^, and with a prefactor of 
= 1.45. The analogy with convection predicts 
that K4 = 27r_fi'i. The small aspect ratio convective 
experiment extrapolated at large aspect ratio gives 
Ki ~ 0.75 X 0.31, which translates into K4 = 1.46. 
This is in very good agreement with the prefactor mea- 
sured by Wendt. 

The second regime predicts torque varying more 
slowly than R^^^. It could therefore only marginally 
explain the second regime observed by Wendt, for 
R > 10'^, in which G ~ R^ "^. However, it could 
explain the regime obtained by [3 [T^] for R ^ 10^, 
in which a continuously varying scaling exponent was 
obtained. This is shown in Fig. [3j where the fit 
to the data of [T5] is compared with the theoretical 
formula (21). The comparison is made using coeffi- 



cients extrapolated from the small aspect ratio con- 
vection experiment: Kq = 20, = 2'kK2 with 
K2 — 0.75 X 0.45. It may happen however that 
this regime 2 does not exist in Taylor-Couette experi- 
ments. Indeed, since the temperature analog is related 
to the velocity, it might be impossible to excite veloc- 
ity fluctuations without exciting pseudo-temperature 



fluctuations. This would mean a direct transition from 
regime 1 (mean flow dominated) to regime 3 (fluctu- 
ation dominated). This possibility is explored in Fig. 
|4] where we show the best fit of the measurements of 
Lewis and Swinney, with formula ( 22 ) . This fit uses 
Kj = 0.50 and Ks = 10^. Notice the big difference 
between these constants and their extrapolation from 
the convective case = 0.145 and Kg = 2x 10''. This 
may reflect the sensitivity to boundary conditions of 
the regime 3. Note however that the fit is excellent 
from R = 10^ up to R = 10^. Below R = 10"^, the 
regime 1 with K2 — 1.46 fits the data very well also. 
As a last check, we have compared this regime 3 with 
the constant fitted for Lewis and Swinney's data, to 
the data of Wendt. The result is shown in Fig. [5] for 
3 different gap i] = 0.68, 0.85, 0.935. The agreement is 
excellent. 
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FIG. 3: Torque vs Reynolds in Taylor-Couette experi- 
ments. The symbols are the data of [TH]. The line is 
the theoretical formula obtained in the hard turbulence 
regime and computed using the analogy with convection 
G = AR^^^/{ln{R/B)f/^. The two constants A and B 
are not fitted to the data, but are analytically computed 
using the analogy with convection. 



F. Velocity fluctuations 

The analogy can also be used to predict the be- 
havior of velocity fluctuations. In [13] . the azimuthal 
turbulent intensity ig — < >/Ue was measured 
at midgap with hot film probes. Above 1 x 10^, a fit 
yields 



ie = O.lOi? 



-0.125 



(23) 



Using the analogy, this intensity is related to the tem- 
perature fiuctuations at mid-gap, in the ultra hard 
turbulent regime (regime 3). The total analog tem- 
perature fiuctuation in fact also includes vertical ve- 
locity fluctuations (see Table 1). In an axisymmet- 
ric turbulence, one could therefore expect that the 
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FIG. 4: Torque vs Reynolds in Taylor-Couette experi- 
ments. The symbols are the data of |13| . The lines are 
the theoretical formula obtained in the soft and ultra-hard 
turbulence regimes and computed using the analogy with 
convection. Soft turbulence eq. (201 (full line); ultra-hard 
turbulence eq. (221 (dotted line) . In the former case, all 
the constants are analytically computed using the analogy. 
In the latter case, we have seek the best adjustment with 
data by adjusting the two constants. 



turbulent intensity measured by Lewis and Swinney 
is proportional to the temperature analog. Recent 
measurements of this quantity at Rayleigh number 
up to Ra = 10^^ have been measured by |2Qj in 
a low aspect ratio Helium experiment. They found 
e'/AT 0.37i^a-°•"^ but this was obtained in a 
regime where the Nusselt number varies like in regime 
2 (velocity fluctuation dominate but NOT tempera- 
ture fluctuation) . Using the analogy, this would trans- 
late into a regime where ~ in clear contra- 
diction with the data of Lewis and Swinney, see Fig. 6. 
This might therefore be another proof of the absence 
of the regime 2 in Taylor-Couette experiment. 

Unfortunately, we are not aware of temperature 
measurements in convective turbulence in the ultra- 
hard regime. In previous analysis of temperature fluc- 
tuations in the atmospheric boundary layer, Deardoff 
and Willis [2] showed that temperature fluctuations 
follow the free convection regime 



Nu 



AT [PrRaNuy/^'' 



(24) 



where the proportionality constant is of the order 1. 
Fig. 6 shows the application of this scaling to the 
data of Lewis and Swinney, where the analogy was 
used to translate torque and Reynolds into Nusselt 
and Rayleigh. The best agreement with the experi- 
mental fit of Lewis and Swinney is obtained for a pref- 
actor 1.8. We can also compare the results with the 
theoretical prediction given by the convective model. 
In this model [23 , the temperature fluctuations in the 



FIG. 5: Torque vs Reynolds in Taylor-Couette experi- 
ments for different gap widths rj = 0.68, r] = 0.85 and 
rj = 0.935. The symbols are the data of [T^. The lines are 
the theoretical formula obtained in the soft and ultra-hard 
turbulence regimes and computed using the analogy with 
convection. Soft turbulence eq. (201 (full line); ultra-hard 
turbulence eq. (221 (dotted line). There is no adjustable 
parameter in this comparison, all the constants being fixed 
either by the analogy with convection, or by the compari- 
son with the data of [13| . 



boundary layer obey: 

Ar2 " "(i?aPr)i/2 l + (ziVu)2 



(25) 



where A„ is the height of the viscous velocity layer. 
The value at the height of the boundary layer is 
obtained for z = Xbl, the size of the bound- 
ary layer, which was shown to vary like \bl ~ 
{RaNu)~'^/^ / ^J\Ti{RaNu). Assuming that the value 
at mid-gap equals this maximal value, we obtain: 



AT 



(26) 



This prediction is shown in Fig. [6] using a fitted pref- 
actors of = 0.16 and Kn = 1. It is in very 
good agreement with the experimental fit of Lewis and 
Swinney. 



IV. CONCLUSION 

In this paper, we have shown how a well-known 
analogy between stratified and rotating shear flows, 
for axi-symmetric perturbations, can be extended into 
the turbulent regime. Assuming predominance of the 
axi-symmetric perturbation in the turbulence dynam- 
ics , we used this analogy in the unstable case (analogy 
between convection and centrifugally unstable Taylor- 
Couette flow) to predict the scaling of the momentum 
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FIG. 6; Azimuthal velocity fluctuations in Taylor-Couette 
flow. The circles are the power-law fits of experimental 
measurements by [T3]. The triangles are the power-law 
fit to the temperature fluctuations (analog of azimuthal 
velocities) in Helium by [20] ■ The line is the prediction 
obtained with the analog of the free-convective regime [2] . 
The dotted line is the theoretical formula predicted by our 
model eq. ( 26 1. 



It would therefore be interesting to see whether the 
analog of the Grossman and Lohse formula, namely: 
G = ci = i?e^/2 + C2i?e^/3 would not fit the data 
equally well than the Eckhardt et al formula. This 
would reduce the number of unknown parameter by 
one. 

It would now be interesting to study in more de- 
tails consequences of the analogy in the stable case 
(i.e. stably stratified flow v.s. centrifugally stable 
flow). There are many observational, numerical and 
experimental results in the case of stably stratified 
flows. However, their counter part in the centrifu- 
stratifled shear flow rotating shear flow 



X 

w 



rdrO, 
2^sin2 

u 

(rv — w cot q 



TABLE I: The detailed analogy between startified and ro- 
tating shear flow. The notations for the stratified case are 
from EH. 



transfer and velocity fluctuations. Our prediction is 
that at low R eyno lds number, the non-dimensional 
torque follows (|20| while at i? > 10", it follows E^. 



The analogy can also be used to discriminate be- 
tween theories about Taylor-Couette turbulent quan- 
tities. For example, we have shown that the "classi- 
cal" Nu ^ Ra}/^ regime, translate into a G ~ Ra^/^ 
in the Taylor-Couette flow (both being unobserved ex- 
perimentally at large Rayleigh or Reynolds number). 

The analogy also sheds new light on the recent 
theory of Eckhardt et al jT7]. It predicts a depen- 
dence: G = cii?e3/2+5«/2 + c2i?e2+3«, where ^ = 
—0.051 is a parameter which has been adjusted to 
a best fit. When translated using this analogy, this 
formula would give in the convective case: Nu = 
CiRa^'^+'^^/'^ + C2i?ai/2+3«/2. This has to be com- 
pared with the theoretical prediction of Grossman and 
Lohse [16 , made using the same theory, which leads 
to Nu = ciRa}/^ + C2Ra}^^- Clearly, there is no 
value of ^ which can reconciliate the two formulae. 



gaily stable rotating case is presently missing. Recent 
experiments by Richard et al |10| performed on flows 
between counter-rotating cylinders could help filling 
this gap. 

Finally, the analogy is of great interest for as- 
trophysical and geophysical applications. In astro- 
physics, for example, many objects are differentially 
rotating, and are characterized by very large Reynolds 
number. These Reynolds numbers cannot be reached 
in laboratory experiments. On the other hand, we 
have at our disposal a natural high Rayleigh (and 
Reynolds number) laboratory of stratified turbulence: 
the atmospheric boundary layer. We believe that we 
could use all the data collected in our atmosphere to 
get great insight about large Reynolds number behav- 
ior of rotating, astrophysical shear flows, using the 
analogy sketched in the present paper. 

We thank Frangois Daviaud for comments on the 
manuscript. 
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